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Abstract 

In a previous paper |l| we performed an analysis of asymptotic structure of perturbation 
theory series for semileptonic r-lepton decays in massless limit. We extend our analysis to 
the Cabibbo suppressed AS = 1 decay modes of the r lepton. In particular we address 
the problem of corrections to theoretical formulas. The properties of the asymptotic 
behavior of the finite order perturbation theory series for the coefficient functions of the 
corrections are studied. 



1 Introduction 



The accuracy of experimental data for r lepton decays makes it feasible now to extract the 
spectral density of the Cabibbo suppressed AS = 1 decay modes through detecting strange 
hadrons and even to pin down the tiny difference with the Cabibbo favored AS = case 
due to the non- vanishing strange quark mass One of the main problems in obtaining 

precise theoretical formulas is the strict control over the convergence of perturbation theory 
(PT) series and the error due to its truncation For the nonstrange decay channels 
(Cabibbo favored) this problem is now an actual problem - the theoretical uncertainty has 
already reached a limiting value existing due to asymptotic nature of the PT series. This 
value is comparable in magnitude with the experimental error. For the Cabibbo suppressed 
modes the experimental errors are still larger than theoretical uncertainties. However, with 
the accuracy of experimental data permanently improving the limiting theoretical precision 
within FOPT is becoming a major problem of theoretical analysis in general and of the 
extraction of the strange quark mass from corrections in particular. 

From the theoretical point of view one of the central quantities of interest for the Cabibbo 
suppressed modes from theoretical point of view is the correction to hadronic spectral density 
arising from the nonvanishing s-quark mass. This makes the description different compared 
to the massless (Cabibbo favored or ud) case. The corrections to the spectral densities 
have been calculated with a high degree of accuracy within perturbation theory in the strong 
coupling constant (e.g. [^). 

In the present note we determine the ultimate theoretical precision reachable for 
corrections within a finite order perturbation theory analysis. We follow closely the lines of 
ref. [jl| and reach our conclusions in a renormalization scheme invariant way. 

The basic observable is the normalized r lepton decay rate into hadrons written in the 
standard form 

-p / It ^ 

Rr = ""i = N,SEWmud\\l + M + \Vus\\l + 5us)) . (1) 

i (r — > Lvv) 
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The leading terms in eq. (0) are the parton model results while the terms 5ud and 5us 
represent the effects of QCD interaction and (in case of nonvanishing quark masses) mass 
effects 10, H, 0, H, H- Vud and Vus are matrix elements of the weak mixing matrix and Sew 
describes the electroweak radiative corrections to the r-decay rate. 

In general, hadronic observables in the r system are related to the two-point correlator 
of hadronic currents with well established and simple analytic properties - this makes the 
comparison of experimental data with theoretical calculations very clean. This feature makes 
T lepton physics an important area of particle phenomenology where theory (QCD) can be 
confronted with experiment to a level of very high precision. 

The correlator (here we concentrate only on strange hadronic current, i.e. the term 
proportional to Vus) has the form 

U^M = 12vr'2 / dxe'^^{Tj^{x)jlm = g^g.n,(g2) + g^,U,{q^) (2) 

with j^{x) = u'~f^{l — 75)5. Ilq{q'^) and Uglq"^) are invariant scalar functions. We work 
within QCD with three light quarks and do not consider corrections due to heavy quarks 



(c-quark) that would enter in higher orders of PT through loop effects [10|. The correlator 
is normalized to unity in the leading parton model approximation with massless quarks. 

The theoretical expression for the QCD part of the decay rate into strange hadrons is 
given by (with Nc\Vus\'^Sew factored out) 

with Rq{s) and Rg{s) being the absorptive parts of the structure functions Hqiq"^) and Ug^q"^). 
The masses of the light quarks {u, d) can be neglected. We study corrections for the 
Cabibbo suppressed decay modes. The representation of the total decay rate in terms of the 
absorptive parts of the structure functions nq(g^) and Tlg{q'^) is convenient from the point 
of view of their analytic properties in the complex g^-plane. The physical decomposition of 
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the correlator reads 

n^.(9) = {q,.q. - g'^M-)nT(g') + g^g.ni(g2) (4) 

where the nr(g^) part contains only spin 1 contributions and 11^ (g^) contains only spin 
contributions. The relation between the two sets of invariant functions IIt,l{(1^) and Uq^glq"^) 
describing the correlator eq. (0) reads 

nT(g^) = ^^#, U,{q') = U,iq^) + ^^. (5) 
—q^ q^ 

In terms of the physical (definite spin) invariant functions eq. reads 

The longitudinal part of the spectral density Rl{s) vanishes if all quarks are assumed to be 
massless. On expanding Ilq{q^) and Ilg{q^) in m1/q^ and keeping only the leading term in 
this expansion one has 

2 

n,(g2) = n(g2) + 3^n^,(g2) (7) 

n,(g2) = -q'Tl{q') + ^m^n„,(g2) (8) 

where n(g^) is the invariant function already known from the mass zero case. The functions 
n^((5^) with = —q^ are computable in perturbation theory in the deep Euclidean region 
^ oo. The results of the PT calculation read 
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Even though the fourth order MS scheme coefficient ks and the coefficients kq2, kg^ are not 
known at present we retain their contributions since we want to dispose on them as free 
parameters for later considerations. The numerical values of the other coefficients in the MS 



scheme are given in the Appendix: these are needed as reference numbers for our transforma- 
tions to other more appropriate schemes. The running couphng and mass are renormahzed 
at the scale /i. The hght quarks u, d are taken to be massless. Eq. (|^) constitutes the 
complete theoretical information necessary for our fixed order perturbation theory analysis 
of ml corrections. The corresponding expressions with an explicit dependence can be 
found by inserting the expansion of the running coupling constant 

TT TT TT Vr Z TT 

+{P,L + ?>(5o(5,L' + hlL' + ^-(3l(3,L^ + (3',L%^f + . . . (10) 

Z 6 TT 

and that of the running mass 

= l + L7o(^) + (iLVo7o + ^^^o^ + ^7i)(^)^ + ... (11) 

Here /3j and 7, are the appropriate coefficients of (3- and 7-functions describing the evolution 
(running) of the coupling and mass and 

^ = M^)- (12) 

The coupling constant in eq. (|T0| , 11) in taken at a genuine normalization point /x. In 
the present note we do not systematically discuss non-perturbative effects stemming from 



standard power corrections |TT[. The standard power corrections arise from nonvanishing 
vacuum expectation values of local operators within the operator product expansion and are 
relatively small. They can be simply accounted for if necessary. They do not mix with the 
corrections. The coefficient functions of the local operators are known in low orders of the 
perturbative expansion and there is no necessity for a thorough analysis of their convergence 
properties at present. 
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2 Natural strange quark masses for internal 

perturbation theory description of mj. corrections 

We restrict our attention only to the new features that appear due to the mass corrections 
when compared to our previous analysis [|l|. The appropriate quantities to consider are 
moments of a spectral density 

M„ = (n + 1) p{x)x''dx . (13) 
Jo 

We often use Mr as a unit of mass which leads to dimensionless variable x = s/M^. Note 
that within finite order perturbation theory the moments eq. ([ISD coincide with the results 



of a contour integration ||I2|, |13[ [Ij, |15| because of analytic properties of the functions In^ z. 
The moments of the hadronic spectral densities are internal characteristics of the hadronic 
decays of the r system and it is instructive to describe these moments in internal variables. 
In the massless limit within FOPT there is only one independent internal variable - the 
effective coupling a{s) which is defined directly on the physical cut through the relation 

p{s) = l + a{s) (14) 

and studied in ref. M. All the constants that may appear due to a particular choice of 



the renormalization scheme are absorbed into the definition of the effective charge e.g. ||16 
[T7| , |18|, |l9l. In our present analysis the effective charge is determined by the massless piece 
n(g^) of the correlator @ (see, eq. (0,11)). Eq. ([T^ fixes the definition of the effective 
charge which is later used as an expansion parameter for the mass corrections. With such a 
definition one retains consistency in the description of the massless approximation for strange 
and non-strange modes. 

The running of the coupling a{s) defined in eq. (0) contains logarithms of s with coeffi- 
cients given by an effective /5-function 

p{a) = = -(3oa^ - (3ia^ - p2a^ - ^ga' + O(a^) (15) 

(J/O 

6 



and reads 

a{s) = a + pola^ + {(3il + /3o'/')a' + (^2/ + ^A/^o/' + /5o^')«^ 

+ (^3/ + 3/3o^2/' + \fill^ + y/Jo'A/^' + /3o'^')«' + ■ ■ ■ (16) 

where a = a{M^), I = ln(M^/s). Note that the expansion of a{s) in eq. ([I6| ) has the same 
form as the one for 0:^(5) /tt but now and jS^ are coefficients of the effective /5-function 
while /3o and (3i are renormahzation group invariants. The effective couphng a can be 
expressed through the MS couphng constant 



a = ^ + 1.64 [-^ j - 10.28 j + (-155.0 + ^3) (^^ j + ... (17) 

with as = as(M^). The massless spectral density reads |]l| 

p{s) = 1 + a + 2.25a2/ + a^(4/ + 5.063/2) + a^(-25.7/ + 22.5/^ + 11.4/3) 

+a^((-409.5 + 4.5A;3)/ - 149.4/^ + 87.75/^ + 25.63/^) + O(a^). (18) 

At any fixed order of perturbation theory the effects of running die out for the high order 
moments (large n in eq. (|13])) improving the convergence of the perturbation theory series. 
With the definition of the charge according to eq. (|1^ all high order corrections vanish 
as n ^ 00 for any fixed order of perturbation theory. However, for the mass corrections 
this is not true anymore since the strange quark mass introduces a new parameter. In 
order to obtain only logarithms of energy in the spectral density one can redefine the MS- 
scheme quark mass appropriately and absorb all remaining constants into the internal mass 
parameter. Because there are two invariant functions that characterize the correlator (H) the 
definitions of the mass parameters may be different for them. This redefinition is nothing but 
the change of subtraction scheme for the given spectral density. Note that the introduction 
of a natural internal coupling parameter such as the effective charge a{s) allows one to extend 
the perturbation theory series needed for the description of relations between observables by 
one more term as compared to the analysis in e.g. the MS-scheme (e.g. ^ 0)- The 
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same reason of obtaining one more term in the perturbation series is behind the definition 
of internal parameters for the quark mass. 

2.1 The ^f-part of the correlator: function Ilmg{q'^) 

First we consider the (yf-part of the correlator which contains spin 1 contributions only. We 
use the massless case as a base for our analysis here. We define a new mass parameter rrig 
for the strange quark to absorb all constants in the mass correction of the spectral density 
into the new mass definition. By definition 

ml{s)pf{s; as) ^ ml{s; a) = ml{Ml)p,{s). (19) 

The new mass parameter = rriglM^) is related to the MS-scheme mass parameter = 
ms(M^) through 

ml = ml(l + 1.67a - S.STa^ - 51. Oa^ + (-1342.5 - 1.67^3 + kg3)a^ + 0(a^)) (20) 




The spectral density Pg{s) for the g-pait of the mass correction contains only logarithms in 
fixed order PT expansion and reads 

Pg{s) = 1 + 2a/ + 0^(8.05/ + 4.25/2) + a3(5.3/ + 38.23/2 + 9.21/3) 

+a^((-45.3 - 2.0k3)l + 67.4/^ + 134.2/^ + 20.14/^) + O(a^). (21) 

With the help of eq. ( pOf ) the numerical value of can be determined in terms of the MS- 
scheme mass squared with an accuracy of about 7% if eq. ( pOD is evaluated up to third 
order. If the series is (perturbatively) inverted, the MS-scheme strange quark mass can be 
extracted with an accuracy of about 4% from the relation 

ml = m2(l - 0.185 + 0.107 + 0.037 + (0.19 + 0.00025A;3 - 0.00015^^3)). (22) 
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This accuracy is sufficient at present for phenomenological applications. 

Given the expression for the spectral density eq. (|21|) the whole analysis of ref. |l|| applies. 
The moments of the spectral density Pg{s) in general behave worse than in the massless case. 
This is understandable if one compares the coefficients of the logarithms in the spectral 
densities (|18|) and (^l]) . The coefficients are larger overall in the mass correction case. The 
basic objects one needs for the construction of observables are moments of the spectral 
density Pg{s) 

Mg{n) = {n + l) f pg{s)s''ds . (23) 

J 

We find 



M,(0) 


= 1 


+ 2a + 16.6a^ + 


137.0a^ + (1378.5 - 2.0A;3)a^ 


M,(l) 


= 1 


+ a + 6.15a2 + 2 


8.67a3 + (141.97- 1.0A;3)a^ 


M,(2) 


= 1 


2 . 
+ -a + 3.63a2 + 
3 


12.31a3 + (35.69 - 0.67A;3)a^ 


M,(3) 


= 1 


+ -a + 2.54a^ + 
2 


6.97a^ + (11.58 -0.5A;3)a^ 


M,(4) 


= 1 


2 . 
+ -a + 1.95a^ + 
5 


4.56a3 + (3.55 - 0.4^3)0^ 



MJlOO) = 1 + a + 0.081a2 + 0.060a3 + (-0.434- 0.0198A;3)a^ (24) 

101 

Note that the unknown coefficient kg-^, which would appear in the fourth order coefficient 
of the moments in the MS-scheme, is absorbed in the definition of the mass m^. Still the 
fourth order coefficient is not known because of its dependence on which enters due to the 
charge redefinition. Indeed, this dependence has its origin in the definition of the mass 
(P!P|). The third coefficient of the effective 7-function 7^3 depends on k^. This dependence 
affects the fourth order coefficient of Pg{s) through the running mass (pT]) 

70 = 1, 7,1 = 4.027, 7,2 = 2.65, 7,3 = -22.65 - ^3 . (25) 

For large n the moments behave better because the infra-red region of integration is sup- 
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pressed. Note that the coefficients of the series in eq. (|2^ ) are saturated with the lowest power 
of logarithm for large n for a given order of perturbation theory, i.e. they are saturated with 
the highest coefficient of the effective /3-function and 7-function. 

It is instructive to compare the results in eq. (0) with the MS-scheme expansions given 

by 



MS I 



M. 



1 + 3.67— + 20.0 

TT 

1 + 2.67— + 6.32 

TT 

1 + 2.33— + 2.70 

TT 

1 + 2.17— + 1.06 

TT 

1 + 2.07— + 0.14 

TT 



+ 110.1 

- 38.98 

- 64.26 

- 73.18 

- 77.46 



vr 

TT 

as_ 

TT 

as_ 

TT 
TT 



+ 



+ 



-256.3 + kg^] 
-1779 + kgs) 
-1865 + fcgs) 
-1863 + fcgs) 
-1851 + kgs) 



TT 
TT 
TT 
TT 
TT 



M^^S(IOO) 



1 + 1.69— + 2.99 

TT 



-87.15 



+ (-1755 + k. 



93 J 



The advantage of the effective scheme against the MS-scheme is apparent starting with 
coefficient. For moments larger than one (ra > 1) in the MS-scheme the series shows already 
asymptotic growth in the third order while in the effective scheme the third coefficient is 
still smaller than the previous one. From its construction it is clear that the convergence 
behavior improves with higher moments in the effective scheme while in MS-scheme high 
moments become even worse than the lower ones from the point of view of the structure 
of the perturbation theory series. The anomalously small third order coefficient of the ffist 
moment in the MS-scheme is the result of an accidental cancellation of the contributions of 
logarithmic and constant terms because of the particular choice of the scheme. The general 
discussion of ref. [|I|] now applies. For our numerical estimates we take a = 0.111 as obtained 
from the corresponding value of the MS-scheme charge. In the effective scheme the PT series 
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for the moments read 



Mg{0) = 1 + 0.222 + 0.204 + 0.187 +(0.21 - 0.0003^3) 

Mg(l) = 1 + 0.111 + 0.076 + 0.039 + (0.022 - 0.00015^3) 

Mg{2) = 1 + 0.074 + 0.045 + 0.017+ (0.0054 -0.00010A;3) 

Mg{3) = 1 + 0.056 + 0.031 + 0.010 + (0.0018 -0.000076^3) 

Mg(4) = 1 + 0.044 + 0.024 + 0.006 + (0.00054 - 0.000061A;3). (26) 

With the choice = 100 |jl[ the fourth order correction is smaller than the third term for 
the moments with n < 5. A formal accuracy of about 0.7% can be obtained if the zero order 
moment is excluded using as an estimate the contribution of the smallest term. With the 
standard Pade estimate for ^3 = 25 the moments become 

Mg{0) = 1 + 0.222 + 0.204 + 0.187 + 0.20 

Mg{l) = 1 + 0.111 + 0.076 + 0.039 + 0.018 

Mg{2) = 1 + 0.074 + 0.045 + 0.017 + 0.003 

Mg{3) = 1 + 0.056 + 0.031 + 0.010 - 0.0001 

Mg(4) = 1 + 0.044 + 0.024 + 0.006 - 0.001. (27) 

If one excludes Mp(0) an accuracy of better than 2% can be obtained. 

There is no value for which makes the forth order corrections of all moments smaller 
than the previous correction. If k^ has a value between —46.5 and 1.2 all moments starting 
from the first moment show no asymptotic growth in fourth order. Such a fine tuning of the 
unknown coefficient k^ seems to be unrealistic. We thus conclude that asymptotic growth 
is unavoidable in fourth order in the g-pait. The ultimate accuracy depends on the value 
of k^ varying between 0.7% and 2% if the zero order moment is excluded. The invariant 
statement about the asymptotic growth is that the system of moments Mg{n) with n = 
included cannot be treated perturbatively at the fourth order of perturbation theory for 
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the given numerical value of the expansion parameter a = 0.111 if one wants to obtain an 
accuracy of the coefficient function in front of correction in Uglq"^) amplitude better than 
15% - 20%. This statement about the ultimate accuracy of the set of moment observables 
attainable in fourth order of perturbation theory is independent of whichever numerical value 
takes. 

The perturbation theory expansions for the system of moments with (1 — s)" weight 

M,(n,0) = (n + 1)^ p,(s)(l - srds = (n + 1)1 ^ ['^^ _ j^y M,{k) (28) 

show a much worse PT behavior. One has 

Mg(l,0) = 1 + 0.333 + 0.332 + 0.336 + (0.397 - 0.00046A;3) 
Mg{2,0) = 1 + 0.407 + 0.429 + 0.461 + (0.569 - 0.00056A;3) 
Mg{3,0) = 1 + 0.463 + 0.509 + 0.572 + (0.728 - 0.00063A;3) . (29) 

The convergence of the series is obviously quite poor. All moments Mg{n,0) contain the 
contribution of Mg{0, 0) = Mg{0) in the sum eq. ( pHD which by itself shows a bad behavior; 
the rest makes it worse. 

To summarize, the general PT structure of the moments for the g-pait in the effective 
scheme with the new mass parameter is very similar to the massless part and there are 
no qualitatively new features found here as compared to the analysis of the massless part in 
ref. 0. 

2.2 The g-part of the correlator: function Ilmq{q'^) 

The g-amplitude Uqlq"^) contains contributions of both spin 1 and spin final states. The 
correction to the g-part is different from the g-pait as concerns its analytic properties in 
the g^-plane - it contains a 1/g^ singularity at the origin - which necessitates a separate 
treatment. The explicit power singularity at the origin of the function Umqiq"^) / q"^ 
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makes the formulation of the moments for the spectral density of the g-part directly on 
the physical cut a bit tricky. Indeed, because of the factor the amplitude for the ml 
correction has no standard dispersion representation. Rather the dispersion representation 
should be written as 

- J JTq~' ^^"^ 

with a measure da{s) which is not differentiable, i.e. da{s) ^ a'{s)ds with some continuous 
(t'(s). However, it can be written in a more familiar form if a different weight is used 

nmg((9') _ r PF{s)ds _ d ^^^^^ ^^^^ 



g2 Jo (s + g2)2 ^g2 

with 



nQ') = L (32) 



r'^ PF(s)ds 
lo s + Q 

and Pf{s) a continuous spectral density. Therefore F{Q^) is the primary function of 



It reads 

Fm = - / ^n^.(g^) ^JdL n^,(g2) L^in^^. (33) 

The result for F{Q^) is 

F(g2) ^L+dL + L'^)a + (15.757L + T.llL^ + lA17L^)a^ + ... (34) 
For the discontinuity across the cut defined by 

Pf{s) = 7^ {F{-s - lO) - F{-s + zO)) (35) 

one obtains 

7 

Pf{s) = 1 + (- + 20a + (1.77+ 14.22/ + 4.25/^)a^+ (36) 
(-207.04 + kq2 + 62.21Z + 54.52/^ + 9.21Z^)a^ + ... 
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where we have aheady substituted the effective couphng a for a^. Eq. ( p6D has the standard 
form of the spectral density for needed for the comparison with the massless and the g cases. 
In order to get rid of the constants the new mass parameter rUg for the amphtude UmqiQ^) 
is defined in analogy with the g-case such that 

ml = mlpF{Ml) (37) 

with TTig = mq{M^) and = m^(M^). Then one explicitly has 

ml = m^(l + ^a+1.77a2 + (-207.044+A;,2)a^ + (-1335.5-2.33A;3-4.92A;,2 + A;g3)a^ + ...) (38) 

Note that the order term contains not only the unknown coefficient kq2 but also the higher 
order coefficient fcgs which makes the contribution of the term completely arbitrary. The 
definition of this new mass in terms of the MS-scheme mass eq. (^B]) has an accuracy of 
about 2% if only second order corrections are used. The PT series is only known to third 
order - already the fourth term contains the unknown coefficient kq2- Expressed through the 
MS-scheme coupling constant as eq. (|38D reads 




With this new mass ml a new spectral density Pq{s) can be defined in analogy to eq. (|T9D 

mlpF{s) = mlpg{s) (40) 

that leads to the expansion 

Pq{s) = l + 2/a + (9.55/ + 4.25/2)0^ + (36.36/ + 44.6/2 + 9.21/3)a3 + 

((-1141 - 2h + 6.75kq2)l + 253.6/2 + 154.96/^ + 20.14/^)a^ + O(a^) . (41) 

Note that the series (^) contains only logarithms and no constants. The coefficients of 
the spectral density Pq{s) from eq. ( |il| ) are close to those of the g-pait mass correction 
Pg{s) eq. (|21| ) up to the second order. The coefficients of the highest powers of logarithms 
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are simply equal while the lower powers are different because of different coefficients in the 
corresponding D-functions eq. (^. The coefficient of the third order term of eq. (^) for Pq{s) 
is seven times larger than the corresponding coefficient from eq. (0) for Pg{s) for the lowest 
power of logarithm which dominates the behavior of the higher moments. The coefficients 
of the spectral density Pq{s) are in general larger than the coefficients of the spectral density 
p{s) in the massless case eq. (|18D . Thus, the behavior of the spectral densities {p,pg,pq) 
allows one to immediately conclude about the convergence of the corresponding moments 
for all three independent contributions up to order in the strange decays. 

The standard moments of Pq{s) as in eq. ([T^j23| ) do not however coincide with the physical 
moments of the g-amplitude. Indeed, the physical g-moments M^^{n) are defined through 
the contour integration in the following way 

where the effective mass is used for normalization. The moments eq. (|^) can be evaluated 
through the function Pq{s) directly. The zero order moment turns out to be 

(0) = -p,(M^2) = -1 . (43) 



The reason for this is that the integration with n = in eq. (|4^ ) picks up exactly those 
contributions that eventually are absorbed into the strange quark mass redefinition. Higher 
order physical moments are related to the standard moments of Pq{s) via 

^rWln>o = n j\q{s)s^-Ms - 1 = Mq{n - 1) - 1 . (44) 

They contain no parton model contribution. The moments Mq{n) are the standard objects 
defined in eq. (|l^) and the whole analysis of ref. |l| is applicable with Pq{s) to be compared 
with Pg{s) and p(s). Numerical values for the moments are 

Mq{0) = l + 2a + 18.1a2 + 180.8a^ + (779.4- 2fc3 + 6.75fcg2)a^ 
Mg(l) = 1 + la + 6.90a2 + 47.39a3 + (-297.3 - /cs + 3.375A;g2)a^ 
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Mg{2) = l + -a + 4.13a2 + 24.08a3 + ( 
Mg(3) = 1 + + 2.92a2 + 15.53a3 + ( 

M„(4) = l + -a + 2.25a2 + 11.28a^ + ( 
5 

Mg(lOO) = 1 + — a + 0.095a2 + 0.37a^ + (-11.25 - 0.020A;3 + 0.067A;q2)a^ (45) 

The use of Pg(s) is a universal way of describing the moments in FOPT and easy for com- 
parison with the massless and the g parts. The physical g- moments M^^{n) are related 
to it by eq. (^) and (Q). They do not, however, change the pattern of PT convergence. 
The moment M^^(l) related to ^^^(0) shows asymptotic growth in the third order already. 
As expected higher moments show a better behavior because the low energy region of the 
integration is suppressed. 

3 Order corrections to r-lepton decay observables 

After introducing the technique for analyzing the moments at order we now apply it to 
an analysis of physical observables. 

3.1 Total decay rate 

The r decay width is given by a specific linear combination of moments. The weight function 
contains the overall factor (1 — s)^ which impairs the convergence of the total decay rate 
observable. The (1 — s)^ factor enhances the infra-red region of integration, i.e. the relative 
magnitude of the contributions of logarithms ln(M^/s) at small energy. The concrete shape 
of the weight function with the weight factor (1 — s)^ is the main source of slow convergence 
of the ml correction to the rate 
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-283.6 - 0.67A;3 + 2.25A;,2)a^ 
-237.13- 0.5^3 + 1.69A;,2)a^ 
-199.699 - 0.4A;3 + 1.35A;g2)a^ 



We retain the different mass definitions for tlie two invariant functions in order to explore 
tlie structure of tlie PT series and to clieck on its convergence. 
Tfie result for tlie mass correction of the total decay rate reads 

Rmr = 6^(Mf(0)-2Mf(l) + Mf(2))-6^(M,(0)-M,(l) + -M,(2)) 

= -6^(2M,(0)-M,(l))-6^(M,(0)-M,(l) + -Af,(2)). (47) 

It is expressed through the effective mass parameters ^ and the moments introduced 
earlier. In the parton model approximation all the moments are normalized to unity which 
makes a glance analysis of eq. {^7\j easy. The convergence pattern of all the moments has 
been obtained already before. Numerically one has 

2 

R^^ = -6^ (l + 3a + 29.21a2 + 314.3a^ + (1856.1 - 3.0^3 + 10.13^52)0^) 

T 

2 

-2^ (1 + 3.67a + 34.84a2 + 337.3a3 + (3745.2 - 3.67A;3)a^) (48) 

The reason for the bad convergence of (^5]) is the contribution of the low moments Mq{0), 
Mq{l) and Mg(0), Mg{l) to the mass correction of the total decay rate. Both series in eq. ( ^8] ) 
converge only marginally calling for the resummation of the series. The total contribution 
is dominated by the g-part which has a three times bigger coefficient in the leading term. 

Inserting the expressions for m^, m^, a in terms of the MS-scheme parameters and 
as into eq. (^) one obtains the standard result 

R^, = -8^ (1 + 5.33^ + 46.0 f + (283.55 + 0.75A;,2) (^Y + ..] • (49) 



TT \ TT / \ Vr 



3.2 The '1+0' Method 

In ref. the numerical value for the strange quark mass in MS-scheme has been extracted 
with the '1+0' method which uses the representation of the total decay rate as a sum of 



(L + T) and L contributions (compare with eq. 

Rr = :^h(l-^Xi(l + 2^] n.,,)(g^) - 2^n,(g^)| ^ , (50) 



,2 \ 2 r / ^2 \ „2 
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n(L+T)(g^) = n,(g2). (51) 

It is assumed that in the MS-scheme the series for the (L + T) part of the correction 
in eq. ( ^0|) converges well The convergence is not impressive in FOPT but the numbers 
given for contour-improved FOPT in ref. show fast convergence. The quantity of interest 
is now 

Bit- = ^/2('l-4^V7l + 24') 3 



2 

Tfl 

6]^^(^r(0) - (2) + 2Mf (3)) 

9 

-6^(3M,(l)-2M,(2)). (52) 



If we look at the moments which are contained in R^'^ ( |5^ it is natural to expect a conver- 
gence behavior better than in the total decay rate because no zero order moments are needed 
to construct R^^- This is an invariant (scheme independent) reason for better convergence: 
the particular combination R^^ receives smaller IR contribution from integration along the 
cut and therefore is better computable in PT. Still the convergence is rather slow. Numerical 
results for the (L + T) part are 

r>L+T\ _ 
^•mr \q— scheme 

2 

- 6^ (l + 1.67a + 12.448a2 + 94.01a^ + (-324.629 - 1.67^3 + 5.625fc„2)a^ + . . .) . (53) 

The convergence persists and the last term is still smaller than the third for the standard 
values of 25 < ^3 < 100 and < kq2 < 160. The total contribution of first three terms is 
0.45 which is a reasonable change of the leading order term due to PT correction. In the 
MS-scheme this becomes 

R^flm-.^.^. = -6^, [1 + 4.0 {^)+ 24.67 [^)\ (-62.77 + k„)[^)\ 

(-3110 + 7.29«:,2 + ^,3) (^)' + ... ) (54) 



TT 



The total contribution of the first two terms amounts to 0.65. The change in the leading 
order prediction is even larger than the total change in the effective g-scheme where one 
more term of PT expansion is available. The only advantage of the '1+0' amplitude from 
PT point of view is the absence of the zero order moment Mg(0) which is the most divergent 



one. Still, the moment M^(l) from eq. (|45|) which contributes to the rate expression eq. ( ^2[) 
is bad enough to prevent a fast convergence. The longitudinal part can be expressed in terms 
of the moments and convergence is rather bad 



TTl Tfl 1 

= -12^(M,(0)-2M,(l) + M,(2))-6^(M,(0)-M,(l) + -M,(2)) (55) 

because again the lowest order moments enter. The leading order parton model term in the 
g-part vanishes. Numerical values are 

4 

2 

- 8^ (O + a + 12.57a2 + 165.20^ + (1635.6 -h + 3.375A;g2)a^) (56) 



i?^^ = - 2^ (l + 3.67a + 34.8a2 + 337.3a^ + (3745.2 - 3.67^3)0' 



M2 

from which one can perceive the reason for the bad PT structure. While the convergence 
of the gf-part is rather standard, the admixture of g-contribution without the leading term 
is enhanced by relative factor of four which makes the sum in eq. (^) completely uninter- 
pretable. 

Note that we use ^ mass parameters for q parts of the correlator as internal mass 
scales for corrections. One can introduce another set of parameters ^ related to the 
definite spin decomposition of the correlator eq. @). As one can see from eq. (1) uit = rUg 
because the H-Tiff) amplitude is proportional to Ilg{q^). The effective mass m\ is obtained 
from th longitudinal part lii^q^). The expression of through the MS-scheme mass is 
reasonable (eq. (pOD) while the corresponding relation for mz, is much wilder. 
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) 



3 



+ 



(57) 



The different structures of PT serieses for the parameters m'^ ans m| can be interpreted 
as a result of the difference of full QCD interaction in spin 1 and spin channels. The 
stronger and earher breakdown of PT behavior in the spin channels can be related to the 
contribution of non-perturbative effects (instantons) which is absent in spin 1 channel. 

To summarize, the convergence of the PT series for the corrections for the most 
natural and precisely measured physical observables is always slow and almost marginal. 
This is because the IR region of integration is numerically important for the given value of 
the coupling constant and the order of PT. Any apparent fast convergence is the result of 
either a specific linear combination of moments or a particular scheme choice. The former 
case is, however, not realized for physical observables of interest measurable in experiment. 
This would imply that a resummation of the series is necessary for a sound interpretation of 
the theoretical formulas for decay rates. 

4 Conclusions 

We have analyzed the asymptotic structure of PT series for corrections. Using the 
standard estimate of the accuracy of an asymptotic series we have found that the theoretical 
precision in the perturbative description of Cabibbo supressed r-lepton decays is already 
limited by the asymptotic growth of the coefficients in fourth order of perturbation theory. 
This is a scheme invariant statement. The accuracy of the perturbative expansion for the 
coefficient functions of the corrections in Cabibbo suppressed channels are 15%-20% 
at best. Therefore the extraction of the numerical value for the strange quark mass from 
the corrections to the r-decay rate into strange hadrons is limited by the precision of 
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the coefficient functions. Better theoretical accuracy can be obtained by using observables 
which contain higher order moments but the experimental accuracy of them is not good 
enough at present. From phenomeno logical point of view the analysis of the corrections 
differ from the massless case. While in the latter the low order moments can be excluded 
by substituting experimantal results for them (from e+e~ annihilation, for instance), the 
coefficient functions of corrections have no immediate physical meaning and cannot be 
traded for in this manner. 

The introduction of two natural mass parameters allows to describe massless, q and g 
parts of the correlator up to order with only two unknown parameters and kq2 instead 
of four in the MS-scheme. The existence of two different mass scales is physically motivated 
by the difference of interaction in spin 1 and spin channels. Still for both independent 
corrections {g and q part) the convergence of low order moments is slow. The contribution 
of IR region is large and these quantities have limited precision when evaluated within PT. 
It can not be improved by a particular change of the scheme or taking particular linear 
combinations. Only those moments converge well where the IR contribution is suppressed. 
The renormalization group improved QCD parameters a(s), ms{s) run too fast to give precise 
PT series for a set of r observables containing low moments with large IR contribution. If 
the running would be slower - coefficients of (3- and 7-functions would be smaller - then 
one could meet the high standards of experimental precision for low moments with FOPT 
theoretical formulas. Observables with higher moments are described well within FOPT but 
experimental accuracy is not yet sufficient for precise comparison. 

Therefore for precise comparison of theory with experiment some procedure of resum- 
mation is required fl^, And, in fact, one can perform resummation only for lowest 

order moments to keep things close to the standard FOPT. Resumation, however, introduces 
additional (to the standard renormalization group freedom [^) arbitrariness, in particular, 
it interferes with nonperturbative power corrections that makes separation of contributions 
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within the OPE not unique. 

This problem is beyond the scope of the present analysis. 
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Appendix 



The input for all calculations are the coefficients of the correlator, the beta and the gamma 
function which are available up to four loop calculations in present. In the case of Ilq{Q^) 
only corrections up to are known because the constant term of the correlator is not yet 



calculable. All coefficients are given in the MS-scheme (see, e.g. p6[). Below ({z) is the 
Riemann (^-function. 



The coefficients of the correlator are ||2^, 

299 



ki 
^92 



24 -9^(2)' 
58057 779^,^, 75^,^, 



432 54 

144 2 ^^^^ 
1967833 vr^ 
5184 36 
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1™5.C(3)-h51!5c(5). 
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108 



Beta function coefficients 



/33 



9 ^ . ^ 3863 

-, di =4, ^2 = • 

4' ^ ' ^ 384 ' 

140599 445,,^, 
+ — C(3). 



4608 



32 



Gamma function coefficients 



7o 



31 



91 

^^ = 24^ 



8885 



73 



2977517 
41472 



+ 



, 72 

32 



|C(3). 



578 



432 
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